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ABSTRACT 



Hydromagnetic dynamo theory provides the prevailing theoretical description for the 
origin of magnetic fields in the universe. Here we consider the problem of kinematic, 
small-scale dynamo action driven by a random, incompressible, non-helical, homoge- 
neous and isotropic flow. In the Kazantsev dynamo model the statistics of the driving 
flow are assumed to be instantaneously correlated in time. Here we compare the results 
of the model with the dynamo properties of a simulated flow that has equivalent spatial 
characteristics as the Kazantsev flow but different temporal statistics. In particular, 
the simulated flow is a solution of the forced Navier-Stokes equations and hence has a 
finite correlation time. We find that the Kazantsev model typically predicts a larger 
magnetic growth rate and a magnetic spectrum that peaks at smaller scales. However, 
we show that by filtering the diffusivity spectrum at small scales it is possible to bring 
the growth rates into agreement and simultaneously align the magnetic spectra. 

Subject headings: dynamo - magnetic fields — magnetohydrodynamics — turbulence 



1. Introduction 



Magnetic fields are thought to be present in a variety of astrophysical objects, including stars, 
planets, accretion discs, galaxies and the interstellar and intergalactic media. In some cases the 
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magnetic field has been directly detected and its properties and consequences have long been known 
(e.g. the Sun's large-scale magnetic field that is responsible for the 22-year sunspot cycle), while the 
active role of magnetic fields in other astrophysical phenomena has only recently been recognised 
(such as its role in angular momentum transport in accretion discs). While in some cases it is 
possible that the magnetic field could be the remnant of a decaying fossil field, in others it is 
believed that the field is sustained against the action of Ohmic dissipation via the operation of a 
hydromagnetic dynamo - a mechanism through which the kinetic energy of the plasma is converted 
into magnetic energy. Understanding the evolution of magnetic fields in the universe remains one 
of the most important unsolved problems in astrophysics. 

Due to the enormous scale of astrophysical objects, in virtually all cases it is believed that 
the flow that leads to dynamo action is turbulent. The generated magnetic field then exhibits 
complex spatial and temporal characteristics on a wide range of scales, a complete understanding 
of which represents a formidable problem. Progress with astrophysical dynamo theory is therefore 
typically made by breaking the problem down into a number of fundamental sub-problems. A 
natural starting point is to consider the initial stages of evolution of a weak seed field. During 
this kinematic phase the field is assumed to be so weak as to have no dynamical effects on the 
turbulence. The Navier-Stokes equations governing the evolution of the flow then decouple from 
the induction equation governing the magnetic evolution, viz. 

— = Vx(uxB) + i|V 2 B, V-B = (1) 

where r] is the magnetic diffusivity. The problem therefore reduces considerably to solving equations 
([PJ, given a flow u(x, t). If the flow is a kinematic dynamo the magnetic field grows exponentially 
and the task is to determine the magnetic growth rate and the preferred lengthscale of the dynamo 
instability. In reality, the exponential growth cannot continue indefinitely. Eventually the gener- 
ated magnetic field will be strong enough to react back on the flow (via the Lorentz force in the 
momentum equation) and the growth will saturate. In this nonlinear phase the aim is to understand 
the dynamics of the saturation process and to determine the amplitude of the resulting field. It is 
important to note a further distinction that is often made: Large scale dynamos are those in which 
the magnetic field has a scale of variation much larger than the characteristic velocity correlation 
length, while in small scale dynamos the magnetic field varies on scales comparable with or smaller 
than the scale of the driving flow. 

In this paper we shall be concerned with the kinematic, small-scale dynamo problem. There 
are essentially three different avenues of research that have developed for addressing this scenario. 
They differ in the way in which the flow is specified. In the first, one considers the amplifica- 
tion of magnetic fie l ds in a laminar flow that has a given analytical form (see the monograph by 



Childress fe Gilbertl (119951 )). The second method is similar, except that the flow is specified nu- 
merically. It is typically either a solution of the Navier-Stokes equations or it is synthesized to have 
certain properties. In the third approach, the flow is assumed to be a random process with pre- 
scribed statistics. In this case the aim is then to determine the corresponding statistical properties 
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of the magnetic field. The work presented in this paper concerns a comparison of the second and 
third approaches. 

The third avenue of research mentioned above (that of considering field a mplification i n tem - 
porally random velocity fields with given statistical properties) was pioneered by lKazantsevi (|1968l ). 
who considered the case of a Gaussian distributed random flow with zero mean that is instanta- 
neously correlated in time and homogeneous and isotropic in space. In this case it follows that the 
flow is completely specified by the two-point, two-time velocity covariance tensor 



(ui(x,t)uj(x + r,t + r)) = Ky(r)5(r), 



(2) 



where (...) denotes an ensemble average and 5(r) is the Dirac delta function. We note that such 
a statistical d e script ion is also often referred to as a Kraichnan statistical ensemble, in reference 



to iKraichnanl (| 19681 ) who studied the evolution of a passive scalar field in a turbulent flow with 
the same statistical properties. The Gaussian property implies that all odd moments vanish and 
all even moments can be expressed in terms of products of t he above s econd order moments. The 
most general form of the isotropic tensor Kjj(r) is (see, e.g.. iBatchelorl ()1953l )) 



Kij(r) = K N (r) ( 8ij ) + n L (r)Af + g(r)e ijk rk, 



(3) 



where 5%j is the unit diagonal tensor, e^-fe is the unit alternating tensor, r = |r| and summation over 
repeated indices is assumed. For an incompressible flow it follows that kn = kl + (r/2)(dnL/dr). 



Assuming that the magnetic field is also homogeneous and isotropic, we have 
(Bi(x, t)5,-(x + r, t)) = M N {r, t) (<% 



^ + M L (r, t) T -^f + K(r, t)e ijk r k , 



(4) 



with the solenoidal constraint implying Mn = Mi+{r/2){dMi J /dr). Interestingly, it is then possible 
to derive a closed system of equations for the evolution of the magnetic correlati on functions [Mr, 
and K) in terms of the corre sponding fu nctions for the velocit y (kl and g) (|Kazantsevl Il968l : 
Vainshtein &: Kichatinovl I1986T ) . Recently, iBoldvrev et al.l (|2005l ) established that the equations 
have a self-adjoint structure, implying that the eigenvalues of the system are real. In the absence 
of kinetic helicity (g = 0) the Kazantsev model predicts zero magnetic helicity (K = 0) and the 
equation for Mi reads 



dM L d 2 M L 
dt dr 2 



+ 



, 4 

K H K 

r 



dM L 
dr 



k" + -k' 
r 



M L 



(5) 



where k(t) = kl(0) — kl{t) + 2rj. Formally, the eigenvalue problem for Ml(t, t) = M(r) exp(Ai) is 
to be solved in r £ [0, oo]. In the non-helical case (the case to be considered here) it can be shown 
that the eigenvalues form a discrete spectrum. The eigenmodes are spatially bound and satisfy 
the boundary conditions M(0) = const, (where for the kinematic dynamo problem the constant is 
arbitrary and is related to the value of t he see d field) and M(r) decays exponentially to zero as 
r — > oo (see, e.g., iMalvshkin fc Boldvrevi (|2007l )). 
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The Kazantsev dynamo model is a very useful tool. It reduces the problem of solving a 
partial differential equation for the evolution of the three-dimensional vector field B(x, t) into 
two ordinary differential equations for the evolution of the scalars Ml(t, t) and K(r,t), given the 
corresponding quantities for the flow. In the former case, computational power limits studies to 
moderately turbulent systems with the Reynolds numbers R e = ul/v (where v is the viscosity) and 
R m being approximately equal and of the order of a few thousand. The validity of extrapolating the 
results to astrophysical situations where the magnetic Prandtl number (P m = R m /R e ) is typically 
either tiny (e.g., in the solar convection zone or the Earth's core) or huge (e.g., in the ISM) and 
the Reynolds numbers are enormous is unknown. By contrast, the Kazantsev equations can be 
solved relatively inexpensively and Reynolds numbers up to 10 12 are permissible (the limit being 
only that of the double prec ision accuracy of the standard computer floating-point format (see 
Malvshkin & Boldvrevl f|2008i )). 



The Kazantsev model is the only available analytical model for dynamo action. It has there- 
fore beco me very popular and over the years the theory has been built upon considerably. For 
example, iKulsrud Anderson! (119921 ) developed a detailed theory for the spectra of magnetic fluc- 
tuat ions driven by an incompressib le velocity field with a Kolmogorov spectrum E(k) ~ A; -5 / 3 (see 
also iKraichnan fc Nagaraianl (|1967l )). The spectr al theory was extended to describe d-dimensional 
flows with arbitrary degrees of compressibility bv lSchekochihin et al.l (|2002i ). The Kazantsev model 



has b een solved for various different forms of the velocity correlator (see, e.g., IZeldovich et al 



(|l990l ) and re ferences therein) and the theory has been broaden e d to a ddress the role of helic - 



ity (see, e.g., 



Bold vrev et al 



Vai nshtein Sz Kichatinovl ()1986l ): iBerger fe Rosnerl (1995); Kim fc Hughes! 

flows with varying degree s of roughn ess (e.g. iBoldyrev fc Cattaned d 



1997) 



([2Q0l|)). 



2004) 



Malvshkin Sz Boldvrev ( 2010 )). nonli near effects ( Kiml 1999) and the consequenc e s on dynamo ac 



tionof finite time correlated flows (e.g. Chandran ( 1997 ) JRogachevskii &; Kleeorin (jl997l ): Schekochihin fc Kulsrud 



Ultimately, as with any simplified model, one must assess the extent to which the solutions of 
the model depend on the assumptions made, in this case regarding the statistics of the velocity. In 
particular, it is important to determine how accurately the Kazantsev dynamo model describes the 
situation in which the velocity is physically realistic, i.e. when the flow satisfies the equations of elec- 
trically conducting fluid dynamics. In this first paper we accept the conditions of incompressibility, 
homogeneity and isotropy, and we also note that the assumption of Gaussianity is not believed to 
be of crucial importance. Indeed, direct numerical simulations have shown that a purely Gaussian 
velocity field amplifies ma gnetic fluctuations in a similar manner to a Navier-Stokes velocity (see 
Tsinober &: Galantil (|2003l )). However, we believe that it is important to note that in addition to 
the magnetic diffusivity, the only inputs to the Kazantsev model are the spatial correlation func- 
tions {kl and g) of the velocity. All of the temporal properties of the flow are encapsulated by the 
5-function in equation ([2]). Indeed, the fact that the flow is instantaneously correlated in time is 
essential for the derivation of the closed system of equations governing the evolution of Ml and K 
(i.e. equation ([5]) in the non-helical case). If the advecting flow is finite-time correlated, as it will 
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be in reality, these equations cannot be obtained in closed form. 

Physically, one might expect that the results of the (5-correlated model would smoothly match 
onto the short-time correlated case if the velocity correlation time is much smaller than the char- 
acteristic time for dynamo action. However, the question remains as to what happens as the 
correlation time of the flow increases. Indeed, it has recently been demonstrated that for the quasi 
two-dimensional dynamo problem, two flows with the same velocity spectrum b ut differing phase 



properties can have very different dynamo properties (jTobias fc Cattaned l2008al ) . The goal of the 
present work is therefore to compare the characteristics of magnetic field generation in the instan- 
taneously correlated Kazantsev flow with the dynamo properties of flows that have a similar spatial 
structure but that are solutions of the forced, incompressible, Navier-Stokes equation and hence 
have finite correlation times. 

Before proceeding, we would like to note that in the remainder of the paper it will sometimes 
prove useful to consider the correlators of the Fourier coefficients of the velocity, rather than the 
spatial correlators ([2]) and ([3]). A further advantage is that the physical meaning of the functions 
kl(t) and g(r) then becomes apparent. We therefore introduce the Fourier transform of the velocity 

u(x, t) = J u(k, t) exp(ik • x)dk. (6) 

In the Kazantsev model, the velocity correlator in k-space then takes the form (assuming incom- 
pressibility) 



{uUKt)u m Qt',t + T)) = {f(A0 U m -^\ +ie lmn k n G(k)\ <5(k' -k)5(r) 



(7) 



where * denotes the complex conjugate. The functions an d q(r) can be ob t ained from 



F(k) and G(k) through three-dimensional Fourier transforms (see iMonin fc Yagloml (|l97ll )). In 
particular, we will use the relation 

. . f°° fsin( kr) —krcos(kr)} . „ 
KL (r) = 8tt / ± — , 3 ~^F{k)dk (8) 

We note that (u(x, t) ■ u(x, t + r)> = 35(t)k l (0) = 25(t) J F(k)dk = 2«J(r) / A-Kk 2 F{k)dk. Similarly 
(u(x,i) • (V x u(x,t + r))) = -28 (t) J 4nk 4 G(k)dk. Thus the functions F(k) and G{k) in equation 
([7]) are related to the kinetic energy and helicity in the Kazantsev model. The equivalent magnetic 
correlator in Fourier space is 

(Bl(k,t)B m (k,t)) = F B {k,t) (s lrn - +ie lmn k n G B (k,t), (9) 

where again M^(r, t) and K(r,t) are related to FB(k,t) and Gs(k,t) through three-dimensional 
Fourier transforms. In particular, EB(k,t) = Airk 2 Fs(k,t) and HB(k,t) = — 87tA; 4 Gb(A:, t) are the 
spectrum functions of magnetic energy and electric current helicity, respectively. 
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2. Formulation of the problem 

We shall first solve the induction equation ([1]) with a prescribed flow that is a statistically 
steady solution of the randomly driven Navier-Stokes equations 

-^ + (u- V)u = -Vp + z^V 2 u + f, V-u = 0. (10) 

Here p is the pressure (whose role is to maintain incompressibility) , v is the viscosity and f is a 
random force, the specific properties of which we detail below. We have neglected the Lorentz force 
since we are considering the kinematic case in which the magnetic field is weak and hence doesn't 
affect the flow. The equations (jl|10p are solved on the triply periodic domain x G [0, 2-7r]. We note 
that although the Kazanstev model is defined on the infinite domain r E [0, oo], it is anticipated 
(and verified below) that the periodic boundary conditions do not play a crucial role if it is ensured 
that the energy containing scales of the flow and the magnetic field are significantly smaller than 2ir. 
It is pointed out that this requirement, together with the restrictions placed by currently available 
computational power, does however severely limit the extent of the inertial range of the simulated 
flow. 

In order to conduct a meaningful comparison with the Kazantsev model the simulated velocity 
must have the required spatial properties, i.e. it must be small-scale (significantly smaller than 
the box size 2ir) and spatially homogeneous and isotropic. For this first study we also restrict our 
attention to the non-helical case. The idea is to impose on the forcing function the properties that 
are required of the velocity. Then, at least in the strongly diffusive case where the solution of 
equation (|10p represents a balance between the diffusive terms and the force, the flow will inherit 
those characteristics. We therefore choose a random, divergence-free, non-helical, homogeneous 
and isotropic force with Fourier coefficients 

f(k,f) = Aa k {k x e k )exp(i0 k ). (11) 

Here A is a (constant) amplitude that is chosen so that the resulting rms velocity fluctuations are 
of order unity, (u 2 ) ~ 1. The k-dependent amplitude, a k) is chosen with the aim of concentrating 
the energy in the velocity fluctuations at scales that are neither too large (since we wish to limit 
the effects of the boundary conditions) nor too small (since the magnetic field will grow on smaller 
scales and we must be able to resolve both and conduct the simulation over a number of eddy 
turnover times at the largest scale). We choose the Gaussian profile 

where k = |k| and we take yo = 7r/4. At each timestep the unit vector e k (i) is rotated about 
one of the three coordinate axes at random by an amount # k (i), with # k (i) and the random phase 
</> k (i) being drawn independently from uniform distributions on [— tt,tt] at each timestep. This 
ensures that the force is short-time correlated. Its correlation time is of the order of the timestep 
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(5t) of the numerical simulations. By decreasing the value of the viscosity we can then vary the 
correlation time of the flow from being approximately equal to that of the force to much larger 
values. This is an important point that we shall return to below. Homogeneity and isotropy of the 
force is maintained by rotating the vector around on the unit sphere. The force is guaranteed 
to be real in physical space (requiring f(— k, t) = f*(k, t)) by initialising e±k appropriately and 
choosing equal seeds for the random number generators for the complex conjugate pairs. It can be 
shown that the spatial correlator of the force takes the form ([3]) with g^ orce {r) = and K^ orce (r) = 
(4/^) (4r 4 /yQ — 28r 2 jy^ + 35) exp(— r 2 /y^), which falls to zero beyond r tt/2 and hence satisfies 
the boundary conditions that the Kazantsev model assumes for the flow. 

For chosen values of v and r] (see below) equations (|l|10p are solved using standard pseu- 
dospectral methods w ith a grid resolution o f 256 3 mesh points (for a detailed description of the 
numerical method see Cattaneo et al. (j2003l )). Initially equation (110p is evolved in isolation until 



the statistically steady state of the flow is reached, which is confirmed by observing the time evolu- 
tion of the energy of the velocity fluctuations. A weak seed magnetic field is then introduced. After 
a short transient evolution the field settles onto the fastest growing eigenfunction and we measure 
the following quantities. 

First, it is necessary to compute the correlation time of the flow, t c say. The spatial correlators 
and spectra defined below represent time averages over a collection of statistically independent 
samples, i.e. samples separated by an interval of the order of the correlation time of the largest 
turbulent eddy. We define r c to be the interval over which the temporal correlator of the velocity 
falls to half of its original value, i.e. the r for which 

(ui(x, t)ui(x, t + r)) t ,v = 0.5(iij(x, t)ui(x, t)) t y. (13) 

Hereafter, angled brackets with subscript t and/or V denote averages over time and/or volume, re- 
spectively. The quantities defined below are computed from data sets that consist of approximately 
50-100 samples separated by an interval of length r c . 

Second, we measure (twice) the exponential growth rate (A) of the magnetic field by fitting 
the function exp(Ai) to (B 2 (x, t))y. We also measure the normalised magnetic spectrum 

EB[k) = 27Tk { <B2(x,t)> v W (W) 

where * denotes the complex conjugate and 9 and (j) are the polar and azimuthal angles of the 
spherical coordinate system onto which the function of k is interpolated in order to construct the 
isotropic function of k = |k|. 

We now wish to compare the growth rate and the magnetic spectrum with that predicted by 
the Kazanstev model. In order to solve the Kazanstev equation ([5]) we require the spatial correlator 
of the velocity k>l(j"), which we compute as follows. First, we calculate the diffusivity spectrum 
from the numerical simulations. We take the trace of equation ([7]) with k = k' and integrate over 
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r, thereby removing the 5-function in time and yielding 

1 f°° 
F(k) = -(u*(k,t) / mQM + T)*-), 



-oo 

« ~«(k,i) / Mi(k,i)dt) w , (15) 

^ Jt=t-nr c 

To obtain the final expression we have replaced the ensemble average with an average over time t 
(assuming ergodicity), replaced the integral over r £ [—00,00] by an integral over a finite interval 
that is sufficiently larger than the velocity correlation time r c in order to ensure convergence (n 
is the smallest integer required for convergence), and integrated over the direction of k. We then 
determine the spatial correlator kl{t) from F(k) by numerically integrating equation ([8]) using a 
modified Simpson's rule0 

We then seek solutions of equation ([5]) in the form M^(r, t) = Ml{t) exp (At). The equations are 
solved numerically using a fourth-order Runge-Kutta method (for details see lMalyshkin fc Boldyrev 
( 200 71)) . The m a gneti c spectrum function Fs{k,t) can be determined from Ml(v, t) by (see 



Monin fc Yagloml fll97l|)) 

1 f°° 

Fs(k, t) = / [r sin(£;r) — r 2 k cos(/cr)]MrYr, t)dr (16) 

4ir z k J 

In the next section we shall compare the Kazantsev magnetic spectrum 2irk 2 Fs{k) and the growth 
rate A with the magnetic spectrum Es{k) and the growth rate that are obtained from the direct 
numerical simulations. 



3. Results 

First, we concentrate on the direct numerical simulations (DNS). We report the results of 
three cases. Case 1 represents a strongly diffusive test case. We take v = 1 in anticipation that the 
velocity will then inherit the properties of the force, especially its short-time correlation. Hence we 
expect that the results of the numerical simulations and the Kazantsev model will be in agreement. 
We choose a small value of the magnetic resistivity, rj = 0.0008, ensuring that the dynamo is 
excited and the magnetic field is resolved. Case 2 represents a more turbulent scenario. We take 
v = rj = 0.007. Thus the nonlinear term in the momentum equation (|10j) is important and the 
turbulent eddies should have a correlation time r c that is of the order of the eddy turnover time 
and is much longer than the correlation time of the force. This particular value of rj {v) is chosen 
to correspond to the case of a marginal dynamo, while ensuring that both the velocity and the 
magnetic field are resolved. Case 3 has v = r] = 0.0025, representing a situation similar to that 



The integration is computed by using a piecewise parabolic interpolation for the non-harmonic functions of k in 
equation (]8]). The result is then multiplied by the sine and cosine functions and integrated analytically. 
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Fig. 1. — The temporal correlation function (ui (x, t)ui (x, t + r))t,v f° r each of the three flows. 
The vertical dotted lines identify the correlation time, i.e. the time at which the correlator has 
decreased to half of its value at r = (see equation (fl3j) ). 



of Case 2 (with the flow being turbulent and the ratio v lr\ = 1) exce 



magnetic growth rate will be larger (see, e.g., ISchekochihin et al.l (|2004l )) 



at that we expect that the 



In all three I 2 ) 1 / 2 w 1.0 and 5t ~ 0.0055. A physical (eddy turnover) time 

unit T p = yo/u rms ~ 0.8 therefore corresponds to a time interval of approximately 145<5t. This is 
to be compared with the correlation time of each the flows, Tj, say, where the numeric subscript 
j = 1, 2, 3 will henceforth denote the case number. Figure Q] illustrates that t\ ~ 105t, r 2 ~ 90<5t and 
T3 ~ 805t. Figure [2] shows the time series of the magnetic energy for each case. The growth rates are 
estimated to be A x w 0.0198±0.0036, A 2 » 0.0257±0.0068 and A 3 w 0.2655±0.0051. The error bars 
are estimated by the bootstrap method, as follows. From the time interval T over which the growth 
rate of the field is measured (i.e. to the right of the vertical dotted lines in Figure [2]) we randomly 
choose I = 1/At sub-intervals of length At, where At is fixed and chosen such that 5t <C At <S T. 
We estimate the field growth rates Aj over each of the sub-intervals, i = 1,2,..., I. The standard 
deviation of the mean rate A = is then estimated as std(X) = std{Xi) / yl . We then 

repeat the procedure a number of times with different At and compute the average {std(X)). It is 
noteworthy that std(X) is not very sensitive to the value of At provided that 5t <C At <C T. The 
error bars are set to ± 1.96(std(A)), which corresponds to the 95% confidence interval (assuming 
asymptotic normality) . 

The solid line in Figure [3] illustrates the magnetic spectra Es{k) obtained from the numerical 
simulations for Case 1. The left-hand panels of Figure H] and Figure [5] show the corresponding 
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Growth rate X=0.01 98±0.0036 Growth rate X=0.0257±0.0068 




50 100 150 200 50 100 150 200 

time t time t 

Growth rate X=0.2655±0.0051 




50 100 150 200 

time t 



Fig. 2. — The magnetic energy as a function of time from the direct numerical simulations: (a) 
Case 1, (b) Case 2, (c) Case 3. The dotted lines mark the beginning of the time interval over which 
the growth rate is measured (see the text). 
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Fig. 3. — A comparison of the magnetic spectrum from the direct numerical simulations (solid 
line) and the Kazanstev model (dashed line) for Case 1. For purposes of reference, the dotted line 
illustrates the diffusivity spectrum 2irk 2 F(k) (n = 12) and the dashed-dotted line illustrates the 
velocity spectrum (defined analogously to the magnetic spectrum, see equation ([HP ). 
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Fig. 4. — Left: As for Figure [3] except for Case 2. Note that the Kazantsev magnetic spectrum 
is shifted to the right by a factor of approximately 2.6. Right: Filtering the diffusivity spectrum 
at kf = 4.85 brings the growth rates into agreement (within the error bounds) and simultaneously 
aligns the magnetic spectra. 



A.(DNS)=0.2655, X(Kaz)=1 .3053 



A.(DNS)=0.2655, A.(Kaz)=0.2661 




Fig. 5. — As for Figure H] except for Case 3. Compared to the direct numerical simulation, the 
Kazantsev spectrum is shifted to the right by a factor of approximately 2.7. Filtering the diffusivity 
spectrum at kf = 5.10 brings the growth rates and magnetic spectra into approximate agreement. 

results for Cases 2 and 3, respectively. We note that the magnetic energy peaks at k\ f» 11, k\ rj 7 
and k\ ~ 12. The dotted lines denote the converged diffusivity spectra, 2irk 2 F(k), calculated by 
using equation (fT5|) . The velocity energy spectra (defined analogously to Es{k)) are shown by the 
dashed-dotted lines. 

We now turn to the solutions of the Kazantsev model. For each of the three cases, using 
relation ([8]) to calculate from the diffusivity spectrum, we solve equation ([5]) for the fastest 

growing eigenmode. We find Ai = 0.0174, A2 = 0.2576 and A3 = 1.3053. The Kazantsev magnetic 
spectra are shown by the dashed lines in Figure [3] (Case 1) and the left-hand panels of Figure H] 
(Case 2) and Figure [5] (Case 3). 

As expected, the agreement between the numerical simulations and Kazantsev results is good 
for Case 1. The growth rates agree within the error bounds and the spectra peak at approximately 
the same wavenumber. However, much less can be said for cases 2 and 3. The growth rates disagree 
by factors of approximately 11 (Case 2) and 5 (Case 3) and the peak in the Kazantsev magnetic 
energy spectra are shifted towards larger wavenumbers (smaller scales) by factors of approximately 
2.6 (Case 2) and 2.7 (Case 3). What is striking however is that the shape of the spectra are similar. 
Indeed, it appears as though the Kazantsev spectrum could be matched to that of the numerical 
simulations by a simple translation. 

In an attempt to determine a possible reason for the apparent translation, we proceed to 
investigate the sensitivity of the Kazantsev dynamo to the energy containing scales of the velocity 
field. We do this by applying a spectral filter to the DNS velocity before we su pply it to the 
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the role of coherent structures on kinematic dynamo action. Here we apply a Heaviside filter to 
the diffusivity spectrum 2irk 2 F(k) such that all modes with k > kf are set to zero. We then assess 
the changes that occur to the Kazantsev dynamo growth rate and the magnetic spectrum as kf is 
decreased. 

It is not surprising that applying the spectral filter at large wavenumbers (i.e. removing the 
energetically insignificant velocity fluctuations) has little effect on the dynamo, however, it is in- 
teresting that the growth rate can be made to pass through the value obtained in the numerical 
simulations by decreasing the filtering scale kf. In particular, we find that filtering at kf = 4.85 
(Case 2) and kf = 5.10 (Case 3) yields the Kazantsev growth rates A2* = 0.0229 (Case 2) and 
^3/ = 0.2661 (Case 3), i.e. brings them into agreement (within the error bounds) with those of 
the DNS. We note that in both cases the optimal filtering scale is very near to the peak energy 
containing scale k\ = 4.5 (Case 2) and k$ = 5. What is even more striking is that for the optimal 
filtering scale the corresponding magnetic spectra are translated to the left so that they almost 
lay on top of those obtained from the DNS. The results are shown in the right-hand panels of 
Figure H] (Case 2) and Figure [5] (Case 3). Thus the filtering procedure simultaneously brings both 
the Kazantsev growth rate and magnetic spectrum into agreement with those obtained from the 
direct numerical simulations. We note that this result is not specific to the cases detailed here. 
We have also checked that the procedure produces similar results when, for example, the flow is 
turbulent and the ratio rj/v is slightly larger than unity. The complementary case of small rj/u 
is also relevant astrophysically but it is inaccessible with the available computational power. We 
have also checked that the results are not specific to our choice of the forcing function, with similar 
results being obtained if the amplitude of the force is replaced by a^/k 2 . 



4. Discussion 

In the Kazantsev dynamo model the statistical properties of the flow are given and one seeks 
the corresponding statistical properties of the generated magnetic field. The Kazantsev flow is 
instantaneously correlated in time, Gaussian distributed with zero mean and homogeneous and 
isotropic in space. In the nonhelical case it is then possible to reduce the dynamo problem to solving 
a single ordinary differential equation for the longitudinal correlation function of the magnetic 
field. This is to be compared with solving the deterministic system, in which a specific flow is 
prescribed and one solves a partial differential equation for the spatial and temporal evolution of 
the magnetic field. The Kazantsev model therefore allows a considerable saving in computational 
effort and permits the exploration of a parameter regime that is much more extreme than would be 
possible otherwise. For these reasons the Kazantsev model has become a popular tool for studying 
astrophysical dynamo action. 

It is necessary however to address the issue of whether astrophysical flows are likely to satisfy 
the model assumptions, i.e. to determine whether a 5-correlated Gaussian random flow is a suitable 
model of the inertial range of astrophysical turbulence. If it is not, we must then determine whether 
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the results of the Kazantsev model depend sensitively on the assumed statistical properties. In this 
paper we have concentrated on the issue of temporal correlations in the flow, as it is believed that 
finite time correlations will naturally arise through the equations of electrically conducting fluid 
dynamics. We have not studied the effects of relaxing the assumptions made regarding Gaussianity, 
homogeneity and isotropy in space and we have also restricted our attention to non-helical flows. 
We have investigated the kinematic dynamo properties of flows with similar spatial structure but 
different degrees of temporal correlations. We have found that while the Kazantsev model accurately 
describes the system when the flow is short time correlated, increasing the correlation time results 
in large differences between the Kazantsev prediction and the evolution of the MHD system, both 
in terms of the dynamo growth rate and the magnetic energy spectrum. In particular, it appears 
that the Kazantsev model overestimates the growth rate and yields a magnetic spectrum that is 
peaked at smaller scales (larger k). 

We believe that we have traced the source of the problem in the Kazantsev model to a mistreat- 
ment of the small, sub-viscous scales of the velocity field. It appears that the model overestimates 
stretching or underestimates magnetic diffusion at small scales. In particular, we have shown that 
the Kazantsev model's magnetic growth rate and the peak in the magnetic spectrum can be simul- 
taneously brought into agreement with the results of the DNS by filtering the diffusivity spectrum 
at scales larger than the scale at which the velocity spectrum peaks. We believe that this is an 
important non-trivial result. 

At this stage we cannot offer a full theoretical explanation of our results. However, we would 
like to suggest the following qualitative physical picture for why the Kazantsev model mistreats 
the small sub-viscous scales and why the filtering procedure works. We note that the Kazantsev 
model assumes that magnetic field is amplified by turbulent eddies that are J-correlated in time 
at all scales. In contrast, in real turbulence the small-scale eddies are not short-time correlated, 
they rotate together with the viscous eddy and have scale-independent correlation times that are 
approximately equal to that of the viscous eddy. As a result, the sub- viscous eddies do not act 
independently to stretch and fold the magnetic field lines as efficiently as the Kazantsev model 
assumes. In the direct numerical simulations presented here the magnetic Prandtl number is chosen 
to be unity and the magnetic energy peaks at the sub-viscous scales (see Figures |4] and [5]) at 
which smooth velocity fluctuations are poorly modelled as a short-time-correlated Gaussian random 
process. By filtering out the incorrectly treated velocity fluctuations at small scales we do indeed 
achieve a better agreement. In addition, it is possible that even better agreement with the Kazantsev 
model could be achieved if the inertial range of our direct numerical simulations could be extended 
(the extremely short inertial range that we currently have is a result of the requirement of limiting 
the effects of the periodic boundary conditions together with the constraints imposed by presently 
available computational power). This issue is certainly worthy of future investigation. 

We would also like to note that a number of other ways of artificially modifying either of the 
two inputs to the Kazantsev model (i.e. r\ and fvc(r)) were tried but none of them were successful in 
the sense of simultaneously aligning both the growth rates and the magnetic spectra. For example, 
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we have found that the growth rates alone can be aligned in a number of ways, for example by 
artificially 'squeezing' F(k) (which is equivalent to rescaling the wavenumber k — > ak, where a is 
a constant), 'squashing' F(k) (i.e. F(k) — > aF(k), where a is a constant), or by simply retaining 
the diffusivity spectrum from the numerical simulations but increasing the magnetic diffusivity r\ 
in the Kazantsev model. However, in all cases the shift in the magnetic spectrum is very slight and 
it does not align with that from the DNS. Filtering the diffusivity spectrum at large scales (small 
k) also didn't work. 

One possible reason why it is difficult to match the Kazantsev predictions with the results of 
the DNS is that, due to numerical constraints, we have a magnetic Prandtl number P m = u/r] 1 
in our simulations. In the small magnetic Prandtl number regime, P m « 1, the Kazantsev model 
may better describe dynamo action since in this regime the magnetic field is generated by a rough 
velocity field that can be treated as a turbulent diffusion process. Indeed, within the inertial range 
of a Kolmogorov flow the correlation time r c ~ I 2 / 1 D\ (where Di ~ vil is the turbulent diffusivity) 
can be estimated as r c ~ / 2 /(/ 1 / 3 /) ~ Z 2 / 3 and is comparable to the turnover time l/vi. In the 
opposite limit, P m » 1, the velocity fluctuations are smooth on the resistive scales and thus 
n{r) = k_l(0) — KL(r) + 2rj ~ — (1/2)k^(0)t 2 + 2r] in equation ([5|). Thus for a given magnetic 
diffusivity, the Kazantzev model is essentially dependent only on the single parameter k^(0). Even 
though in this case the sub-viscous eddies are not short-time correlated, and hence by the above 
arguments the Kazantsev model is not expected to work perfectly, essentially any change to the 
diffusivity spectrum that adjusts k^(0) to the right value in the sense of matching the magnetic 
growth rate must also simultaneously align the magnetic spectra. Unfortunately, numerical tests of 
the small and large magnetic Prandtl number regimes must await further increase in computational 
power. 
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